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^ ■ Abstract 

'^ 
p i' We study stochastic inflation in the presence of a dynamical gravitational con- 

Q . stant. We describe the Arnowitt-Deser-Misner formalism for Jordan-Brans-Dicke 



theory of gravity with an inflaton field. The inflaton and dilaton scalar fields can 
be separated into coarse-grained background fields and quantum fluctuations. We 
compute the amplitude of the perturbations generated by those quantum fluctuations 
in JBD theory with an arbitrary potential for the inflaton field. The effect of the 
quantum fluctuations on the background fields is equivalent to a Brownian motion of 
^ ! the scalar fields, which can be described with the use of a Fokker-Planck diffusion 

equation. The probability to find a given value of the fields in the comoving frame 
can be written as a Gaussian distribution centered on their classical trajectory, with 
decreasing dispersion along both field directions. We also calculate the condition for 
the Universe to enter a self-regenerating inflationary phase. The probability distri- 
bution in the physical frame, which takes into account the expansion of the proper 
volume of the inflationary domains, will be concentrated at the Planck boundary and 
will move along it towards large values of the fields. 



^ E-mail: bellido@slacvm.slac.stanford.edu 



Introduction 



The early inflationary scenarios could not solve but only postpone the problem of initial conditions 
for the Big Bang. The first models [|ll assumed that the universe started in a very hot state that 
supercooled in a metastable vacuum, which then decayed to the true vacuum through a first 
order phase transition or just rolled down through a second order phase transition. Linde's 
chaotic inflation [^ opened the possibility of arbitrary initial conditions for arbitrary effective 
potentials for the infiaton field. It was soon realized that the most natural initial conditions 
for inflation were close to Planck scale, where quantum fluctuations of the space-time metric 
become important. 

With the inclusion of the fields' quantum fluctuations, it was understood [^, Q that the 
whole Universe might be in eternal self-reproduction of inflationary domains, while our own 
observable universe could be just a late subproduct of one of such domains. The Brownian 
motion of the infiaton field under the effect of its own quantum fluctuations can be described 
with the help of diffusion equations, in the context of the so-called stochastic inflation formalism 
0, ^. This formalism describes the interplay between quantum fluctuations and coarse-grained 
background fields, and allows one to study the very large scale structure of the Universe 0-0, 
much beyond our observable universe. It is also a consistent framework for the description of the 
origin of density perturbations during inflation, that may have given rise to the recently observed 
temperature fluctuations in the cosmic background radiation |T^ and, furthermore, to the density 
inhomogeneities for galaxy formation. 

Stochastic inflation stresses the fact that inflation is a random process, in which quantum 
fluctuations of the fields act as stochastic forces on their own background values as a result of 
the expansion in de Sitter space. The universe is divided into causally independent inflationary 
domains, in which the fields acquire different values. The global description of this stochastic 
process uses diffusion equations for the probability distribution to find a given value of the fields 
in a given inflationary domain. The probabilistic description of stochastic inflation opened the 
possibihty of a connection with quantum cosmology [0, ||, |[ , in the sense that it may provide an 
interpretation of the wave function of the universe as the probability amplitude for an ensemble 
of inflationary domains, instead as for the whole Universe. 

It is generally assumed in stochastic inflation that the gravitational dynamics of the whole 
Universe is correctly described by general relativity. However, it seems a strong assumption to 
extrapolate the description of the gravitational phenomena at our local and low energy scales 
to the very large scales beyond our observable universe. In fact, it is believed that the theory 
of general relativity is just a low energy effective theory of the gravitational interaction at the 
quantum level. So far the only consistent but by no means definite, since we lack the experimental 
observations needed to confirm it, theory of quantum gravity is string theory [|1T|. String theory 



contains in its massless gravitational sector a dilaton scalar field as well as the graviton. The 
low energy effective theory from strings has the form of a scalar-tensor theory, with non-trivial 



couplings of the dilaton to matter [|^, |13|. Therefore, it is expected that the description of 



gravitational phenomena close to Planck scale should also include this extra scalar field. 



The simplest scalar-tensor theory is Jordan-Brans-Dicke (JBD) theory |T^ with a constant 
u) parameter, where the Brans-Dicke scalar field is in this case the string dilaton, which acts like 
a dynamical gravitational constant. Jordan-Brans-Dicke theory of gravity was introduced in the 



inflationary universe scenario as extended inflation [|T5], in an attempt to solve the graceful-exit 
problem of old inflation. However, only small values of the BD parameter were compatible with 
observations of the cosmic background radiation, in contradiction with the large values coming 
from post-Newtonian experiments. A possible solution could be a variable uj parameter [16|, see 



however |T^ . Extended Inflation has also been considered in the context of string effective actions 



Tsj . A different proposal is to consider JBD theory together with chaotic inflation [ |T9(| . 



In this paper we study the stochastic inflation formalism in the context of JBD theory of 
gravity with an arbitrary chaotic potential for the inflaton. The presence of a scalar fleld in the 
gravitational sector of the theory complicates the picture. There are new effects associated with 
the quantum fluctuations of both the inflaton and dilaton flelds; the former deflnes the effective 
dynamics of inflation and the latter the fluctuations of the gravitational constant [^ . The main 



idea of stochastic inflation is to solve the equations for the inhomogeneous flelds in de Sitter space 
by separating both the gravitational and scalar flelds in short-distance quantum fluctuations, 
which oscillate on scales smaller than the Hubble radius, and long-distance fluctuations which 
are treated as classical flelds. It is the special character of de Sitter space and its horizon which 
allows this distinction. As the inflationary domains expand, short wavelengths cross the horizon 
and communicate with the long wavelength classical background flelds, acting as a stochastic 
noise term. Therefore we can separate the evolution of the long wavelength fleld as a two-step 
process, flrst a random kick from short-distance quantum fluctuations (diffusion) followed by 
classical evolution (drift). This process can be understood as a Brownian motion of the flelds 
and thus can be described with the Langevin and Fokker-Planck equations. 

In Section ^ we describe the Arnowitt-Deser-Misner (ADM) [^ formalism for JBD theory 
with an inflaton scalar fleld, see also [p2|. In this Section we follow closely the work of Salopek 



and Bond [^ and generalize it to JBD. We flnd the energy and momentum constraints for the 
inhomogeneous flelds, the gravitational and the scalar flelds' equations. These non-linear equa- 
tions are extremely difficult to solve and one usually expands in small perturbations with respect 
to a homogeneous background. Fortunately, in the context of stochastic infiation, one is allowed 
to separate the short-distance behavior from the long-wavelength background fields, which can 
then be described by a coarse-grained homogeneous field. The background field equations are 
then obtained by neglecting large-scale gradients. There is a natural scale in de Sitter space given 
by the Hubble radius. Since we are interested in structures larger than the horizon scale it is 
natural to expand in powers of k/{aH) where k is the comoving wave number of the perturbation. 
This is done in Section ^, where we arrive at a consistent set of background field equations. In 
Section ^ we solve these equations in the synchronous gauge and in the slow-roll approximation, 
where we assume that the scalar fields evolve slower than the expansion of the universe. This 
simplifies significantly the equations and we obtain approximate solutions for the background 
fields for generic chaotic potentials of the type V{a) = ^cr'^"'- In Section | we describe the for- 
malism of small perturbations on a homogeneous background, using linear perturbation theory 



of the metric and scalar fields in the longitudinal gauge p3. We find the amplitude of pertur 



bations in the slow-roll approximation for the Einstein frame, in terms of both the dilaton and 



inflaton fields. Those perturbations have a quantum mechanical origin in the short-wavelength 
fluctuations of the scalar flelds. We compute the amplitude of those quantum fluctuations in 
a de Sitter universe and flnd the Gibbons-Hawking temperature [^ for both the inflaton and 
dilaton in certain well deflned limits. We then obtain an expression for the amplitude of density 
perturbations in the cosmic background radiation that should be compared with the one coming 
from general relativity. 

In Section ^ we describe the stochastic inflation formalism P, ^ P] for Jordan-Brans-Dicke 
inflation, see also [^ . The Brownian motion of the scalar flelds in de Sitter space can be written 
as a Langevin equation for the coarse-grained flelds with an effective white noise generated by 
quantum fluctuations. The associated Fokker-Planck equation for the probability distribution of 
flnding a given value of the scalar flelds in a given point of space-time can be derived from the 
Langevin equation and written in a form which is manifestly time-reparametrization invariant. 
The proper way to flnd stationary solutions is to solve the diffusion equation subject to certain 
well deflned boundary conditions 0. We show that the probability distribution far from the 
Planck boundary behaves like a Gaussian with flelds centered at their classical trajectories and we 
calculate its dispersion coefficients. The probability distribution in the physical frame takes into 
account the exponential growth of each inflationary domain. The condition for self-reproduction 
of the universe is then calculated. Those inflationary domains with flelds inside the region of self- 
reproduction will dominate the proper volume of the universe and will tend to diffuse towards 
the Planck boundary. 

In Section |^ we study the stationary solutions of the Fokker-Planck equation in the physical 
frame by reexpressing the evolution equation as a Schrodinger equation in two dimensions with an 
effective potential. We flnd that the distribution will very quickly settle at the Planck boundary. 
However, since the Planck boundary is a line it will slide along it, unless we impose extra boundary 



conditions |20|. In Section N we present our conclusions. 



2 ADM formalism for JBD Cosmology 



In the Arnowitt-Deser-Misner formalism |^ the four dimensional metric g^y is parametrized by 
the three- metric hij and the lapse and shift functions N and A^*, which describe the evolution of 
time-like hypersurfaces p6|. 



goo = -N'^ + h'miNj, goi = g-Q = Ni, gtj = Kj. (1) 



The action for the Jordan-Brans-Dicke theory plus an inflaton scalar field with potential V{a) 
in the ADM formalism has the form |2^ 



S = / d'^x^—g 
= fd^xNy/h 

271 



,. '-R-lid<j^r-lidar-Via) 

iOTT 2 2 



(2) 



{WY-a\ia\' ~V{(J 



where Mp((/)) = — 0^ is the dilaton-dependent Planck mass and 11''^'= are the scalar- fields' mo- 
menta 

, (3) 

n'^ = -(a-iVV|,). 

Vertical bars denote three-space-covariant derivatives with connections derived from /ijj, ^^^i? 
is the three-space curvature associated with the metric hij, and Kij is the extrinsic curvature 
three-tensor 



Kij = ^{Ni\j + Njii-hij), 



(4) 



where a dot denotes differentiation with respect to the time coordinate. The traceless part of a 
tensor is denoted by an overbar. In particular, 



K 



V 



i'i Q *i ' 



K = Kl 



(5) 



The trace A' is a generalization of the Hubble parameter, as will be shown below. Variation of 
the action with respect to A^ and N"^ yields the energy and momentum constraint equations 






MjJ 



KL - Ik\, + - ((t)\,Kl + \(t>\,K] + ^^m%i + ^^^ {TlU\i + nV|,) = 0. 



"i\i 



(6) 



M2(0) 



Variation with respect to hij gives the dynamical gravitational field equations, which can be 
separated into the trace and traceless parts 



3.-, 



All 



K - N^K, = - Nl: + N (^-c(3)i? + -K,K- + -K^ + ^^^^^ 



T 



2 S S AA 



cP\2 



(7) 






(8) 



Variation with respect to and a gives the scalar-fields' equations of motion 



1 

iV 



(fi'^ - iv^nj) - irn^ - ^iV|,0i* - 0|f + ^^(K - N'K\,) 



Ml{<P) 



AiKpN 



'-^^R + Ki.K'^ + -K^ 



dV 



The energy density on a constant-time hypersurface is 



0. 



1 r 



(n*)2 + 



^N 



+ 



{WY + a\,a\' +V{a), 



(9) 



and the stress three-tensor 



T- 



«j 



'\i^\i + ^N^li + ^ij 



{li^f 



1 

+ 2 



:n'^)' - o\,a\^ 



V{a) 



(10) 
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It is extremely difficult to solve these highly nonlinear coupled equations in a cosmological 
scenario without making some approximations. The usual approach is to assume homogeneity of 
the fields to give a background solution and then linearize the equations to study deviations from 



spatial uniformity. The smallness of cosmic microwave background anisotropics |jTO[ gives some 
justification for this perturbative approach at least in our local part of the Universe. However, 
there is no reason to believe it will be valid on much larger scales. In fact, the stochastic 
approach to infiation suggests that the Universe is extremely inhomogeneous on very large scales 
0. Fortunately, in this framework one can coarse-grain over a horizon distance and separate 
the short- from the long-distance behavior of the fields, where the former communicates with 
the latter through stochastic forces. The equations for the long-wavelength background fields are 
obtained by neglecting large-scale gradients, leading to a self-consistent set of equations, as we 
will discuss in the next section. 



3 Spatial gradient expansion 



It is reasonable to expand in spatial gradients whenever the forces arising from time variations of 
the fields are much larger than forces from spatial gradients. In linear perturbation theory one 
solves the perturbation equations for evolution outside of the horizon: a typical time scale is the 
Hubble time H~^, which is assumed to exceed the gradient scale a/fc, where k is the comoving 
wave number of the perturbation. Since we are interested in structures on scales larger than the 
horizon, it is reasonable to expand in k/{aH). In particular, for infiation this is an appropriate 
parameter of expansion since spatial gradients become exponentially negligible after a few e-folds 
of expansion beyond horizon crossing, k = aH. 

It is therefore useful to split the fields and a into coarse-grained long-wavelength background 
fields (f){t,x^) and residual short- wavelength fiuctuating fields 6(j){t,x^). There is a preferred 
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timelike hypersurface within the stochastic inflation approach in which the sphtting can be made 
consistently, but the definition of the background field will depend on the choice of hypersurface, 
i.e. the smoothing is not gauge invariant. For stochastic inflation the natural smoothing scale 
is the comoving Hubble length {aH)~^ and the natural hypersurfaces are those on which aH is 
constant. In that case a fundamental difference between between and 50 is that the short- 
wavelength components are essentially uncorrelated at different times, while long-wavelength 
components are deterministically correlated through the equations of motion. 

In order to solve the equations for the background fields, we will have to make suitable 
approximations. The idea is to expand in the spatial gradients of and to treat the terms that 
depend on the fluctuating fields as stochastic forces describing the connection between short- 
and long-wavelength components. In this Section we will neglect the stochastic forces due to 
quantum fluctuations of the scalar fields and will derive the approximate equation of motion for 
the background fields. We retain only those terms that are at most first order in spatial gradients, 
neglecting such terms as 0,^ , 0|i0'*, ^^^R, ^^-"-Ri, and T/. 

We will also choose the simplifying gauge A^* = 0. The evolution equation (|^) for the traceless 

- - - • d 

part of the extrinsic curvature is then K^j = NKK^- — '^K'^j4>- Using NK = ——Invh from 

(1), we find the solution K^ ex 0~^ h~^^'^, where h is the determinant of hij. During inflation 
h~^''^ = a~^, with a the overall expansion factor, therefore Kj decays extremely rapidly and can 
be set to zero in the approximate equations. The most general form of the three-metric with 
vanishing Kj is 

hij = a^{t,x'') 7ij(x^), a{t,x^) = exp[a(t,a;'')], (12) 

where the time-dependent conformal factor is interpreted as a space-dependent expansion factor. 
The time- independent three-metric 7jj, of unit determinant, describes the three-geometry of the 
conformally transformed space. Since a(t, x^) is interpreted as a scale factor, we can substitute 
the trace K of the extrinsic curvature for the Hubble parameter 

Furthermore, since the Brans-Dicke parameter u is bounded by post-Newtonian experiments 



28| and primordial nucleosynthesis |^9[ to be very large, u > 500, we will use the approximation 
a; 3> 1 in equations (|§-|^). The energy and momentum constraint equations (^ can now be 
written as 

together with the evolution equation (|^ 



where T = 3 [HU^f + ^{Wf - Via)). 

In general, if is a function of the scalar fields and time, H{t,x^) = H{(f){t,x^),a{t,x^),t). 
From the momentum constraint (|Iip we find that the scalar-fields' momenta must obey 



n^ = ^ ^p^'^) 



An 
r2i 



'dH\ H 



4n \ da 



(16) 



Contrary to the situation in Einstein theory |2^, H has an explicit time dependence. Comparing 
equation (p!5| ) with the time derivative of H{(j),a,t), 



^(dH\ =u^(dH\ +u^(9JL] ^^(dH\ 

(p Mp2(0) V^ ^ ^ ' ^ N\dt 



(17) 



1 / f)fT\ Q 

we find — — — = —HIr'. It is only in the slow-roll approximation, when we neglect 11''^ 
versus H, that the Hubble parameter becomes time- independent. 

On the other hand, the scalar fields' equations (^) can be written to first order in spatial 
gradients as 

1 -^ ^ 3M2(0) „ 
—U^ + SHW^ = ^^-^H^, 

N 2n (^g) 

1 • dV 

N da 

This set of equations are still too complicated to solve for arbitrary potentials V{a). In the next 
section we will find solutions to them in the slow-roll approximation. 



4 Jordan— Brans— Dicke Inflationary Cosmology 



We are now interested in the classical behavior of the long- wavelength quasi-homogeneous back- 
ground fields, and later will study the effect of the stochastic perturbations on these background 
fields. Furthermore, we will chose the synchronous gauge N = 1, together with N^ = 0, and 



solve the homogeneous part of the inflationary equations of motion (|T8]), (|T^ in the slow- roll 
approximation. 



r2 rr2 _ SttV (a) 



H' = HL = ^, (19) 

^^ 3Mp(0)' ^ ^ 



together with <^ Hcf) <^ H'^ip. The field equations then look like |]19 

-if. 



M^{<P)dH 



a 



271 d(j) ijj 
Ml{<p) dH _ V\a) (20) 

47r da ~ 3/7 ' 
d = i7 

Let us now study the behavior of the approximate infiationary solutions for a generic chaotic 

2n^ 




potential, V{a) = ^a^", 



For these theories, and (y9 = J^ a move along a circumference of constant radius in the 
plane {(p,(p)- We can parametrize the classical trajectory by polar coordinates {(f){t),{p(t)) = 
{r sin 6 {t),r cos 6 (t)), with constant radius r, and angular velocity given by 

/ , \ 1/2 n^ 

^^^^=U^ (i) ' t^^^^^W]""'- (22) 



For n = 1 we find a constant angular velocity and the solutions can be written as |T^ 
(j)[t) = r sm I Bo 

T 
o{t) = —j= COS I Qo + 




T^'^'^K^^mr (23) 



<^) = «o I \:^ n 1 ~ ^ ' ^o^ < ^ < ^^' 

with the usual power-law behavior of extended infiation. For n > 2, the angular velocity decreases 
and the classical solutions are more complicated, although still power-law. These solutions are 
actually attractors of the complete equations of motion, for all n 



In the chaotic infiation scenario, the most natural initial conditions for infiation are set at 



the Planck boundary, V{ap) ~ Mp(0p), beyond which a classical space-time has no meaning 
and the energy gradient of the inhomogeneities produced during infiation becomes greater than 
the potential energy density, thus preventing infiation itself. The initial conditions for infiation 
are thus defined at the curve (h'i = ^\ -^(J^. On the other hand, infiation will end when the 
kinetic energy density of the scalar fields becomes comparable with the potential energy density, 

102 + 1^2 ^ y^^^ ^^ ^^ ^ i^Mp(0e) = ^0e. 



In the absence of any potential for 0, the dilaton remains approximately constant after 
inflation, and therefore the Planck mass today is given by its value at the end of inflation, 

Mp ~ y — 0e, while the total amount of inflation is approximately — ~ — . 

V ^ Op \(t)p) 

Now that we have the background solutions we would like to describe the formalism of small 
perturbations on the homogeneous background, to account for the quantum fluctuations in de 
Sitter space. 



5 Perturbations on a homogeneous background 



In order to compute the perturbations on a homogeneous long- wavelength background, we use 
linear perturbation theory of the metric and scalar fields in the longitudinal gauge and apply it to 
our case of Brans-Dicke theory with an inflaton field. We follow the notation of WM , and define 



0(t,X^) = 0(t) + (50(t, xO 

ait.x') = a(t) +5a(t,x') 

. . (24) 

N{t,x')= l + <l>{t,x'), N'{t,x') = 

h,,{t,x')= [l-2$(t,a;0]a'(t) S,, . 

Note that we are calling (j) to the homogeneous part, to avoid carrying unnecessary subindices, 
and 50 to the inhomogeneous perturbation. The extrinsic curvature (^) and the three-curvature 
scalar take simple expressions. 
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ir = -3i7(l - $) + 3$, Kij = 0, (%=-$|f, (25) 

(Jj 

while the scalar fields' momenta are fl*^ = 0-t-(50 — 0$. The homogeneous part of the momentum 
constraint (|^) is trivial, but that of the energy constraint is not. 



^-M^b^+r +^'<^))-2^^- p«' 



The inhomogeneous parts of the energy and momentum constraints have complicated expressions, 
\$J' - 3H'^<^ - 3/7$ = 3H- ( 2$ + — - -i I + 3-$ 



a2 I' 



+ cx6a - 0^$ - 0-2$ - 0^— - 0-2— + -T—Sa - 2V- 



(27) 



Mp2(0) V (j) da 
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The homogeneous part of the gravitational equation (J^) is 



2H + ?>H' 



^^ ^1 12 1 -2 



M2(0) V2 



+ o^ - y{(^) + 2^T- 



while the inhomogeneous part can be written as 



(29) 



$ + AH^ + 2H^ + 3i/^$ = i/- 2$ + — - ^ + -$ + 



47r 



M2(0) 
Putting together ( p6l) and (p9|) , we find 






da 



(30) 



H= --H' 
a 



An 



Ml{4>) 



/)2 + (t2) +AH-. 



(31) 



Furthermore, the scalar field equations can be separated into homogeneous and inhomogeneous 
parts 



(f) + 3H(f)= —H' 

dV 
a + "iEd = - ^—. 
da 



6a + 3H6a 



1 
1 



a; 



5(T,J* - 4(T<i> 



,2 |i 



^12a2 1^ 
da"^ da 



•) 



(32) 



(33) 



Subtracting (|30D from (p7|) and using equations (pSj), (|31|), and (p2D, one finally finds a long 
wavelength solution very similar to that of reference 



$ + 



A 1- 



arft 



(34) 



where A is a constant that can be calculated during infiation, using (pSj) and (^) in the slow-roll 
approximation, as 

a f (j)6(j) + a6a\ 



A 



J adt 



9 + a^ 



(35) 



In the Einstein frame g^j^^ = (p'^g^u, the density contrast of perturbations that enter the horizon 
during the matter era is related to both $ and 6(f) by [|4| | 






'-H 



+ a6a 



9 + a^ 



(36) 



^This expression was first obtained in [Q in the Einstein frame. See also 
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The energy density perturbations we observe in the cosmic background radiation could have 
originated during inflation in our model @, as quantum fluctuations of both scalar fields that first 
left the horizon during inflation and later reentered during the radiation or matter dominated eras. 
The amplitude of energy density perturbations (^6]) associated with a given wavelength should 
be evaluated during inflation at the time in which that wavelength first crossed the horizon. 



5.1 Quantum Fluctuations of Scalar Fields 

We will now compute the amphtude 6(f) and 6a of the quantum fluctuations of the scalar fields 
in de Sitter space, whose wavelengths are stretched beyond the horizon and act on the quasi- 
homogeneous background fields like an stochastic force. In order to calculate this effect we 
coarse-grain over a horizon distance and separate the scalar fields into long- wavelength classical 
background fields (f){x) and a{x) plus short-wavelength quantum fluctuations 6(j){x) and 6a{x) 
with physical momenta k/a > H, 



6(f){x,t) = d k 6{k — saH) akUk{x) + alul{x) 



J, 



6a{x,t) = d k' 6{k' — eaH) bk'Vk'{x) + bl,vl,{x) 



(37) 



where e is an arbitrary parameter that shifts the scale for coarse-graining p. Physical results 
turn out to be independent of the choice of e. The quantum fluctuations are assumed to satisfy 
the following commutation relations 

K, 4,] = [h, bl] = 6%k - k') , [afc, bU = . (38) 



The approximate solutions to the perturbation equations (|33D in de Sitter space with V{(7) 
^m'^a'^ are given by P2 | 

Jk-x u^ nog 



UkiX 



^V^H^:\kv), fi' = - + 



rs^ 



(27r)3/2 2 v'"/^> v-'/y , t^ Au~ A 



'■-w- (l^f^^^^'f*'")- "^-M^-j 



(39) 



\IX] 



where r] = —(aH) ^ is the conformal time, and Hi^L{x) = —J-^ (l + ^) exp(i 

The amplitude of the quantum fluctuations of and a can then be computed as 

,1/2 H 



(A'Kk"\Uk 



3|„ |2^ _ 

~27r' 

./3| |2\l/2 H 



(40) 



6a = [A-nk \vk'\ , - o ' 

zvr 



which coincides with the Gibbons-Hawking temperature p5|. It should be noted that we have 



made two approximations when deriving this formula, first that /i and v are both 3/2, see 
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and second that H is constant, when in general this is not true. We should have corrections to 



but they are expected to be small |23 



Substituting the amphtude of quantum fluctuations of the scalar fields ( ^OD into the energy 
density perturbations (|36D and using the equations of motion (pO]), we find for arbitrary V{a) 



1 + 



5p _ 2AH{a,(t)) V{a) 

J ~ y M2(0) V'{a) [l + {^/a 



a 



(41) 



A^A 



where Nx stands for the number of e-folds before the end of inflation associated with the horizon 
crossing of a particular wavelength. For perturbations of the size of the present horizon, we must 
compute the last expression at A^a ~ 65 0. Note that in the large u limit, 4> '^ a during the last 
stages of inflation, which ensures the approximate equivalence of the Einstein and Jordan frame. 
We then recover the usual expression M 



6p^ 2A H{a,(f)) Vja) 
7 ^ y M2(0) V'{a) 



(42) 



Nx 



where Mp is now 0-dependent. For theories with potentials of the type 2~^ > it behaves like 



Quj /2cjA\^^V"+i 



^- — I- 



(43) 



Nx 



In the case of the theory Ao"'^/4, the density perturbation (^31) takes the usual V^ dependence. 
However, in the case m?a'^ /2, the perturbation can be written as 



5p_ 
P 



2e- 



5V37r 



1 + 3lu 1 



m 



M„ 



(44) 



Therefore, we note that the larger is the Planck mass at the end of inflation in a given region of 
the Universe, the smaller is the density perturbation in this region. For a detailed discussion of 
its cosmological implications see ref. [^ . 



6 Stochastic Inflation 



We can now put together the results of previous sections and study the Brownian motion of the 
background fields under the stochastic impulses of their own quantum fluctuations. We could 
write the evolution of the coarse-grained fields in the form of Langevin equations 



Act = n'^A^At + AS^ , 
Aa = HNAt + ASa , 



(45) 
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where AS^^^ are the stochastic noise terms acting on the long-wavelength background. They are 
Gaussian distributed with zero mean and variance 




(AS^M AS^^it')) = 5i, ^P^ H{<Pk)NAt 5{t - t'), (46) 



which corresponds to the Brownian motion of the background scalar fields and a with step 
(^). It is then straightforward to derive the corresponding Fokker-Planck equation, see e.g. the 



appendix in [^ 



^^^-^^^^cc^^')"^^-: 



dPc_d_ ( MldH 

'dt~ aaV47a^^'^' ' M"^ da , 



- ^ S^'-^^-^^^l^^lf™')"^^. 



27r dcj) 87r2 dcj) ^^ ' > ) d<\) da 

where we have chosen the Stratonovich version of stochastic processes. Pc{a, 0; t) represents the 
probability distribution of finding a given field configuration [a, (p) at a given point (of size H~^) 

dP 



in CO moving space. Equation ( |47| ) can be interpreted as the continuity equation — — h V ■ J = 0, 



associated with the conservation of probability. The first terms of each current correspond to the 
classical drift forces ( POj ) associated with the scalar fields, while the second terms correspond to 
the quantum diffusion due to short-wavelength fluctuations (EO). 



It can be seen by inspection that there is an exact time reparametrization invariant solution 
given by 

P.(a,0;t)~(OT^)-V^exp(^^) , (48) 



which is proportional to the square of the Hartle-Hawking wave function of the Universe |3^ , for 
an arbitrary potential V{a) and a variable gravitational constant. Unfortunately, this function 
is not normalizable for potentials that vanish at their minimum. But the real problem is that 
the inflationary regime ends much before the flelds reach the minimum of the potential, and thus 
equation (^71) is not valid there. Therefore, there are no stationary solutions for P^. in realistic 
inflationary models. The proper way to flnd stationary solutions during the inflationary phase 
would be to consider the distribution Pp and solve the corresponding diffusion equation subject 
to certain boundary conditions [^]. This was done explicitly in ref. |2^; we will return to this 
question in the next section. 

Instead of looking for stationary solutions of equation ( ^Tj ) we will try to understand the 
general behavior of the probability distribution P(.{a,(j)\t) in the synchronous gauge iV = 1, for 
generic potentials V{a) = ^cr'^"'- Pc then satisfles the equation 

3/2 



dPr d I 2uj\ 







dt da \2iu \ 3n J da Svr^ \ 3n J da 



1/2 a, 1 /o, ,\\3/2 ^ \- (49) 



dfh \ I, J \ ^n. I dfh Stt^ \ ?,n. I dfh V 
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In order to follow the evolution of the scalar fields under classical drift forces and quantum 
diffusion processes, it will be useful to introduce some redefinitions, 






(50) 



We will follow closely the work of Goncharov and Linde 0, generalized to JBD. We define the 
new functions A and B as 



3n /27r\2 9/i 



^' = -2-x{-) ^V"^'^*' ^*- X\./ 
and write the evolution equation for the probability distribution. 



{a,(P), B = h\aA), (51) 



dPe d 



du da 
d_ 



-A,{a, <P)P,{a, 0) + B'/\a, 0)^ {b'^^cx, 0)P,(a, , 
-A^cr, 0)P,(a, 0) + B'/\a, 0)^ {B'/\a, <p)P,{a, . 



(52) 



We will now consider the case in which V{a) <^ Mp(0), where quantum diffusion is much 
smaller than classical drift, 



g^B{a,(j)) < \A^{a,(j))\, ^B{a,<P) < |A^(a,0)|, 



(53) 



and the classical solutions to the equations of motion satisfy adu) = A^{ac,(pc) and (pc{u) = 
A^{ac,4>c) respectively. Under a further change of variables, x = a — (rc{s), y = (p — (f>c{s) and 



s = u, the evolution equation (|5^ can be written as 
dPr (dA„, 



ds 



^(^c,0c)j ^[yPc{x,y)] + B{a,,(P,)^, 



(54) 



where A„{a, 0) - A„{ac, (pc) ^ a;— -^ + y^r-^ and A^{a, 0) - A^{ac, 0c) ^ V^tt' + ^^T^- The 
cross-terms have been neglected since they contribute with vanishing correlators. 

Under the initial condition P{x,y; 0) = 6{x — Xo)S{y — yo), the dispersion coefficients can be 
calculated as M 



a 






" Lo Al 3r22 y^Ti) 

A\ ?,n\2TTj M02 (7,2/ 3r2 V27ry 61 



o vc , 



(55) 
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The probability distribution P^ then looks like a Gaussian with decreasing dispersions along both 
field directions, centered on the classical field trajectory in the plane {a, 0), 

P.(a,0;t)~exp|-^^^^^^ 2A|(^) ' ^''^ 

This distribution describes the diffusion of the scalar fields in comoving space after a given 
initial condition, and shows how the fields follow the classical trajectory with variable dispersion. 
However, in this description we are not taking into account the expansion of physical space. 



6.1 Self-reproduction of the inflationary universe 

If all the infiationary domains contained fields that followed their classical trajectory with a given 
dispersion, all would have ended infiation by now. However, it was realized in ref. |^ that those 
few domains that jump opposite to the classical trajectory contribute with a larger physical space 
and therefore dominate the physical volume of the universe. Those domains will split into smaller 
domains, some with lower values of the scalar fields, some with higher. As a consequence of the 
diffusion process, there will always be domains which are still inflating, and this corresponds to 
what is known as the self-reproduction of the inflationary universe. 

We are thus interested not in the probability distribution in the comoving frame Pc, but 
on the distribution Pp{a, 0; t) over the physical volume, which takes into account the different 
quasi-exponential growth of the proper volume in different parts of the universe, 

Pp{a, <P- 1) ^ P,{a, <P- 1) e3^*(^-^^) (57) 

where Pc{cr,(j);t) is given in (^). We can now study the onset of the self-reproduction of the 
inflationary universe. The bifurcation line is deflned by the values of the scalar flelds for which 
the maximum of the probability distribution Pp ( |57D starts increasing, or equivalently, where the 
quantum diffusion of the flelds dominates its classical motion. Its expression is best found in 
polar coordinates, in which the classical motion takes the form, 

^~ (/? ' z u ^ ^ ^' (58) 

r^ = 0^ + ifi'^ — > f = 0, 

while the quantum diffusion steps can be written as 



z 





(1 + .=) 



2vr0 ^ ^ (59) 

rSr = [{(j)6(j)f + {Lp5(p)] = -— r. 

It is clear that along the radial direction quantum diffusion will always dominate classical motion, 
since f = 0. However, along the angular direction the condition that quantum diffusion dominates 
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classical motion during the time interval At = H ensures the self-reproduction of the universe 



only in the range [19 



3Mp^( 



1 + 



< Via) < M'U). 



(60) 



Depending on the value of tu, this could be a very wide range of values for the scalar fields. Once 
an inflationary domain enters this regime it will create more domains with growing values of the 
scalar fields. In particular, the inflaton field will tend to diffuse towards larger and larger values, 
which favor a larger volume of physical space. A reasonable question would be to ask if there is 
a stationary solution to the upward drive of the inflaton due to its quantum diffusion. This will 
be examined in the next section. 



7 Stationary Solutions 



The best way to analyze the possibility of a stationary regime to the upward diffusion of the 
inflaton is to look for stationary solutions of the Fokker-Planck equation for the probability 
distribution in the physical frame. 



dt 



+ 



d_ 
d 



47r da " ^ 87r2 da ^ ' 

Ml{<p)dH H^'^ d 



-Pv + 



(h'/'p,) 



(61) 



+ 3HP„. 



27r d(j)''' ' 87r2 

Due to the term 3HPp there are no simple stationary solutions. However, with the use of the 
identity D'^g — 2D{Df g) = e-^D'^{e~^g) — [{Df)'^ + D'^f]g, we can write the stationary probability 
distribution for an arbitrary potential as 

P,(a,0;t) ex e^*i7(a,0)-^/2exp f^^||^) ^(^,0), (62) 

where the new function \E'(o", 0) satisfies a two-dimensional Schodinger-like equation. 



'i73/2^\ ^ ^ (h'^^j-) ^ - Vi^^ </')^ = 8^'^^. 



(63) 



with an effective potential 

167r^ 
V{a,<P)= ^- 

64:71* 



A-n"^ ( 

E~\aA)V\af ^ —E-\a,4>)[ 



H-%a,(t>)Via) + —H(a, 



4 V{a) 



Via) 



(64) 



247r^if((T, 



3uj uj 

For general theories of the type V{a) = ^cr^", the Schrodinger potential 

2 ^5 



can be written as 



(l^)^'V,.. 



~) \ 4A" 



a 



27r\2 3n2(n + 2) 



iX 



a 



2-n 



+ 



27ry /3n\2 03 /27r\2 /9n\ a" 27r /187™\ a'' 



(65) 
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Even in this case it is extremely difficult to find exact analytical solutions to the Schrodinger 




equation (^), due to the non-linear coupling of the two scalar fields. In ref. ^^ we performed 
a numerical simulation that gave us insight into the behavior of the distribution. We will try 
to give here a general idea by analyzing a simple generic case, n = 2. By transforming into 
polar variables (|58|) , we notice that the Universe will enter the regime of self-reproduction for 
z < (\uj'^/12iT'^y^^ <ti 1 and therefore we can approximate the diffusion equation by 

^iJ + l^L^iJ = Sn'E(^) ^. (66) 

Xu z \Xuj/ 

The same approximation {Xuj^ <^ 1) ensures that the second term, corresponding to the diffusion 
in the angular direction, will dominate over the ffist. Therefore, it is expected that due to 
the potential, the distribution will become concentrated close to the Planck boundary, where 
z ^ Zp = vA uj/4:7i. Once there, the distribution will slide in the radial direction along the 
Planck boundary, see ref. |2^. Under the change of variables r = 4/s^, the diffusion equation 
can then be approximated by 

where E" includes a term coming from the 2;-diffusion. Equation (|67D does not have a stationary 
solution. It is well know from quantum mechanics that a potential proportional to — 1/s^ is too 
strong and produces singular solutions at s = (r = oo). The probability distribution will move 
forever towards large r, i.e. it is non-stationary. This result is generic, it is true for all potentials 
of the type Ao"^". Those probability distributions were named "runaway solutions" in ref. [pO|] . 

Runaway solutions are a special feature of Jordan-Brans-Dicke, since the evolution in general 
relativity is equivalent to the motion along the z direction, which is stationary. The dynamics 
of the gravitational constant introduces the r-motion along the Planck boundary, that is non- 
stationary unless we impose further boundary conditions. For instance, it was shown numerically 



in ref. pO[ that a sharp increase in the inflaton potential could effectively simulate a boundary 



condition in a, and therefore produce global stationary solutions. 



8 Conclusions 



In this paper we have described the mathematical formalism for stochastic inflation in the context 
of Jordan-Brans-Dicke theory of gravity. If string theory is the correct description for the grav- 
itational interaction at the quantum level, it is important to know the behavior of the massless 
gravitational sector of strings in the low energy effective theory, below the Planck scale. Such 
a sector contains a dilaton scalar field as well as the graviton, and therefore it is worthwhile 
studying the cosmology of scalar-tensor theories of gravity. The dilaton acts like a dynamical 
gravitational constant, the Brans-Dicke field. Even though it is argued that the present exper- 
imental observations put severe constraints on such a scalar component of gravity, we cannot 
disregard its effects close to the Planck scale, where quantum fluctuations of the fields become 
important. 



Since chaotic inflation, it was realized that the inflationary regime went all the way up to 
the Planck scale. The quantum fluctuations of the fields then became important and their effect 
on their own background fields was described with the use of diffusion equations. Stochastic 
inflation is the mathematical formalism that describes the Brownian motion of those coarse- 
grained background fields during inflation. It is also a very powerful tool for the description 
of the very large scale structure of the Universe, much beyond our observable universe. Its 
probabilistic description opened the possibility of a connection with quantum cosmology 0. In 
order to incorporate the effect of the string dilaton in the description of the inflationary universe 
close to the Planck scale, we have studied the stochastic inflation formalism in the context of the 
simplest scalar-tensor theory, Jordan-Brans-Dicke theory of gravity. 

For that purpose we described the Arnowitt-Deser-Misner formalism for JBD theory with 
an inflaton scalar field, in order to include in a consistent way the non-linear evolution and 
backreaction of the scalar fields and the metric. We then separated the fields into long- wavelength 
background fields and short-wavelength quantum fluctuations, thanks to the natural scale of the 
de Sitter horizon in inflationary domains, and obtained a self-consistent set of equations for the 
background fields by expanding in spatial gradients. Solutions to these equations are easily found 
in the slow-roll approximation, describing the classical evolution of the scalar fields. 

Quantum fluctuations of the scalar fields during inflation give rise to adiabatic energy den- 
sity perturbations that could be responsible for the temperature fluctuations recently observed 
by COBE in the cosmic background radiation. We describe a consistent way of computing the 
amplitude of the perturbations on a homogeneous background in terms of the quantum fluctu- 
ations of both dilaton and inflaton scalar fields. For theories of the type Acr^ we find the usual 
scale-invariant vA dependence. However, for theories with a mass term for the inflaton, the 
amplitude of adiabatic density perturbations is proportional to the ratio of such a fundamental 
scale (assumed given) and the effective value of the Planck scale at the end of inflation. 

We then studied the new stochastic phenomena that appear when the dilaton is included in 
the diffusion equations for the probability distribution of finding a certain value of the scalar 
fields both in the comoving and the physical frame. While the probability distribution in the 
comoving frame followed the classical trajectories of the scalar fields, with some time- dependent 
dispersion, the description in the physical frame takes into account the quasi-exponential growth 
of the inflationary domains. Those few domains that jump opposite to the classical trajectory 
inflate more and thus contribute with a larger physical space. There is a bifurcation point (in fact 
a line in JBD inflation) at which the contribution to the proper volume of the Universe of a given 
inflationary domain dominates the tendency of the fields to follow their classical trajectories, and 
the Universe enters a regime of self-reproduction. In that regime some domains end inflation and 
enter the radiation and matter dominated eras, while others continue inflating and producing 
new inflationary domains, overwhelmingly dominating the total proper volume of the Universe. 
In that sense, our own observable universe might be an insignificant subproduct of one of those 
domains that ended inflation [^j. We compute the bifurcation line in the case of JBD stochastic 
inflation and find that the range of values of the scalar fields in the self-reproduction regime 
depends on the value of the Brans-Dicke parameter. 
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The main question to be addressed is the existence of stationary solutions for the probabihty 
distribution in the physical frame. With the introduction of the dilaton field, the diffusion space 
is now bidimensional and the self-reproduction of the Universe occurs on an infinite wedge in this 
space, see eq. (|60D, bounded between the bifurcation line and the Planck boundary although open 
for arbitrary large values of the scalar fields. We have found, see also ref. [^] for a more detailed 
description based on numerical simulations, that in general we will obtain runaway solutions: 
non- stationary probability distributions concentrated close to the Planck boundary but moving 
with different speeds towards large values of the scalar fields, in particular to an infinitely large 
Planck mass. Perhaps this is an indication that we need to introduce a dynamical cutoff in the 



inflaton potential. In ref. [BO] we discussed the cosmological consequences of this result. 



In the case of general relativity, the energy eigenvalue of the stationary distributions was 
related to the fractal dimension of the global inflationary universe [^. In our case, the existence 



of runaway solutions indicate a very complicated fractal structure of the universe close to the 
Planck scale. Further work should be done in this direction. 
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